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Gamma-ray bursts as dark energy-matter probes in the context of 
the generalized Chaplygin gas model. 
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ABSTRACT 



In this paper we consider the use of Gamma Ray Bursts (GRBs) as distance markers to study 
the unification of dark energy and dark matter in the context of the so-called Generalized 
Chaplygin Gas (GCG) model. We consider that the GRB luminosity may be estimated from 
its variability, time-lag, and also use the so-called Ghirlanda relation. We evaluate the im- 
provements one may expect once more GRBs and their redshift become available. We show 
that although GRBs allow for extending the Hubble diagram to higher redshifts, its use as a 
dark energy probe is limited when compared to SNe la. We find that the information from 
GRBs can provide some bounds on the amount of dark matter and dark energy independently 
of the equation of state. This is particularly evident for XCDM-type models, which are, for 
low-redshifts (z < 2), degenerate with the GCG. 
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1 INTRODUCTION 



The GCG model jKamenshchik, Moschella & Pasquier 2001] 
|Bento, Bertolami & Sen 2002bT "is' an interesting alternative to 
more conventional approaches for explaining the observed accel- 
erated expansion of the Universe such as a cosmological constant 
(see e.g. lBento & Bertolami 1999l|Bento, Bertolami & Silva 2001 
or quintessence <Ratra & Peeb?eH988al fl988b1 
Wetterich 19881 ICaldwell, Dave & Steinhardt" l998 

Ferreira & Joyce 1998| |Zlatev, Wang & Steinhardt 1999| 

Binetruy 1999T IRim 19991 lUzan 19991 IChiba 19991 
Amendola 19991 lAlbrecht & Skordis 20001 |Fujii 2000) 

Bertolami & Martins 20001 |Sen, Sen & Sethi 2001 1 

en & Sen 200T1 |Bento, Bertolami & Santos 2002) . It is worth 
remarking that quintessence is related to the idea that the cos- 
mological term could evolve iBronstein 1933 Bertolami 1986a 
IBertolami 198rjbl IQzer & Taha 19871 and with attempts to solve 
the cosmological constant problem. 

In the GCG approach one considers an exotic equation of state 
to describe the behaviour of the background fluid: 



< a < 1 is considered. Within the framework of Friedmann- 
Robertson- Walker cosmology, this equation of state leads, after be- 
ing inserted into the relativistic energy conservation equation, to an 
evolution of the energy density as iB ento et al. 2 002b l 
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where A and a are positive constants. The case a = 1 corresponds 
to the Chaplygin gas. In most phenomenological studies the range 
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where a is the scale-factor of the Universe and B a positive in- 
tegration constant. From this result, one can understand a striking 
property of the GCG: at early times the energy density behaves as 
matter while at late times it behaves like a cosmological constant. 
This behaviour suggests the interpretation of the GCG model as an 
entangled mixture of dark matter and dark energy. 

This model has several attractive features. From a theo- 
retical point of view, the pure Chaplygin model (a = 1) 
equation of state can be obtained from the Nambu-Goto ac- 
tion for d-branes moving in a (d + 2)-dimensional spacetime 
in the light cone parameterization ( Bordem ann & Hoppe 1993| . It 
is also the only fluid which admits a supersymmetric general- 
ization jJackiw & Polychron akos 2000 1, and also appeared in the 
study of the stabilization of branes in bulks with a black hole 
geometry l Kamenshchik, Moschella & Pasquier 2000 1. The Chap- 
lygin gas may be viewed as a quintessence field with a suit- 
able potential l Kamenshchik et al. 2001 1, or as an effect arising 
from the embedding of a (3 + 1) -dimensional brane in a (4 + 
1) -dimensional bulk l Bilic, Tupper &"Viollier 2002) . The gener- 
alized Chaplygin gas also has a connection with brane theories 
<Bento et al. 2002bl . The model can be yet viewed as the sim- 
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plest model within the family of tachyon cosmological models 
jFrolov, Kof man & Starobinsky 2002 1. 

The GCG model has also been successfully 
confronted with different classes of phenomenolog- 
ical tests: high precision Cosmic Microwave Back- 
ground Radiation data l Bento, Bertolami & Sen 2003a 
Bento, Bertolam i & Sen 2003b| |B entoH3ert o?amHH$eT 2003c 
Caturan & Finelli 2003 Amendola et al. 2003 1, supernova data 
jFabris, Goncalves & Souza 2002a [ | Dev, Alcaniz & Jain 2003] 



^^^^^^^^^^^^^^^roch ella 2003 | 



Makler, Oliveira & Waga 2003| | Alcaniz, Jain & Dev 2003] 

B^rtolamiTt^i^004l lB^ntoeTaT~200^r and gravitational 
lensing tSilva & Bertolami 20031 |Dev, Jain & Alcaniz 2004) . 
More recently, it has been shown using the latest supernova 
data jTonry et. al. 2003] IBarris et al. 20041 IRiess et al. 20041 , 
that the GCG model is degenerate with a dark energy model 
with a phantom-like equation of state ( Bertolami et al. 2004 
IBento et al. 2005t . Furthermore, it can be shown that this does 
not involve any violation of the dominant energy condition 
and hence does not lead to the big rip singularity in future 
IBertolami et al. 20041 . It is a feature of GCG model, that it 
can mimic a phantom-like equation of state, but without any 
kind of pathologies as asymptotically the GCG approaches 
to a well-behaved de-Sitter universe. Structure formation in 
the context of the Chaplygin gas and the GCG was origi- 
nally examined in I Bento et al. (2002b) , |Bilic et al. (2002) and 
|Fabris, Goncalves & S ouza (2002FJ. The results of the various 
phenomenological tests on the GCG model are summarized in 
IBertolami (2004] and |Bertolami (2005) , 

Subsequently, concerns about such an unified model were 
raised in the context of structure formation. Indeed, it has been 
pointed out that one should expect unphysical oscillations or 
even an exponential blow-up in the matter power spectrum 
(San dviket al~2 003 l. given the behaviour of the sound velocity 
through the GCG. Although, at early times, the GCG behaves like 
dark matter and its sound velocity is vanishingly small, at later 
times the GCG starts behaving like dark energy with a substantial 
negative pressure yielding a large sound velocity which, in turn, 
can produce oscillations or blow-up in the power spectrum. This is 
a common feature of any unified approach when the dark matter 
and the dark energy components of the fluid are not clearly identi- 
fied. These components are, of course, interacting, as they make 
part of the same fluid. It can be shown however that the GCG 
is a unique mixture of interacting dark matter and a cosmologi- 
cal constant-like dark energy, once one excludes the possibility of 
phantom-type dark energy (Bento, Bertolami & Sen 2004 1. It can 
be shown that due to the interaction between the components, there 
is a flow of energy from dark matter to dark energy. This energy 
transfer is not significant until the recent past, resulting in a negli- 
gible contribution at the time of gravitational collapse (z c ~ 10). 
Subsequently, at about z ~ 2), the interaction starts to grow yield- 
ing a large energy transfer from dark matter to dark energy, which 
leads to the dominance of the latter at present. Actually, it is shown 
that the epoch of dark energy dominance occurs when dark matter 
perturbations start deviating from its linear behaviour and that the 
Newtonian equations for small scale perturbations for dark matter 
do not involve any mode-dependent term. Thus, neither oscillations 
nor blow-up in the power spectrum do develop. 

In this paper we study yet another cosmological test and its 
possible use to study the GCG. Schaefer (2002 1 suggested that 
Gamma Ray Bursts (GRBs) may be used to extend the Hubble dia- 
gram to redshifts as high as z ~ 5. For 'ordinary' dark energy, such 



high redshifts are not very interesting since at those epochs the Uni- 
verse is dominated by dark matter, and thus it is less sensitive to the 
nature of dark energy. 

For the GCG however, the GRB test might be relevant since it 
unifies dark energy and dark matter into one single fluid. Therefore, 
within the framework of the GCG model, the dark matter domina- 
tion period actually depends of the nature of the dark energy com- 
ponent that kicks in at later times, and one can expect that the study 
of the matter dominated era will bring some insight on some prop- 
erties of GCG models. 

This paper is organized as follows; in section 2 we ex- 
plain our method of using the time-lag/luminosity and variabil- 
ity/luminosity correlations to constrain cosmological models. In 
section 3 we present and comment the results obtain from this 
method. In section 4 we consider a more precise correlation found 
by Ghirlanda, Ghisellini & Lazzati (2004a l, and study its conse- 
quences. In section 5 we consider whether it is possible to use the 
extended redshift range of GRBs to break the degeneracy between 
the GCG and the XCDM (cold dark matter plus a dark energy com- 
ponent with a constant equation of state). Finally in section 6 we 
discuss our results and present our conclusions. 



2 METHOD. 
2.1 Overview. 

The starting point of our study is the proposed correlation be- 
tween time lags in GRBs spectra and the isotropic equiva- 
lent luminosity I Norris, Marani & Bonnell 2000 1, and the correla- 
tion between GRB variability and isotropic equivalent luminosity 
(Reichart et al. 200H . The time lag, denoted by Ti ag , measures the 
time offset between high and low energy GRB photons that arrive 
on Earth. The variability, V, is easily defined in qualitative terms as 
a measurement of the "spikiness" or complexity of the GRB light 
curve. The isotropic equivalent luminosity is the inferred luminos- 
ity of a GRB if all its energy is radiated isotropically. That is, if P 
is the peak flux of a burst in units of photons cm -2 s _1 between 
observer frame energies Ei and E u , the isotropic equivalent peak 
photon luminosity of the burst in erg s _1 between source frame 
energies 300 keV and 2000 keV is given by 



U ao = 47ir 2 (z)P 



J™°EN[E/(l + z)]dE 
J^N(E)dE 



(3) 



where N(E) is the observer frame spectral shape, usually param- 
eterized by a Band function (Band et al. 1993 1, and r(z) is the co- 
moving distance to a burst at redshift z. 

The possible use of this relation to expand the Hubble diagram 
to higher redshifts was first discussed in Schaefer (2002 1. One lim- 
itation of the employed method is the cosmological distances of 
GRBs, which affect the ability of performing a proper calibration 
of their distance independently of the background cosmology. One 
has to either fix a cosmological model and find a calibration that 
depends on the assumed cosmological model, or fit the data to both 
calibration and cosmological parameters. This degrades precision 
since there are more free parameters for the same number of data 
points. 

One way around this was proposed by Takahashi et al. (2003 1. 
Let us assume that one measures the luminosity distance up to z = 
Zmax, with say z max — 1.5. This is, for instance, very likely to be 
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possible with the SNAP 1 mission. This means that one would have 
an estimate of the absolute isotropic luminosity that is independent 
of the calibration and of the cosmological model. One may then use 
these estimates to calibrate the (r; ag , Li SO ) and (V, Li so ) relations 
without assuming a background cosmological model. 

The major strength of GRBs as cosmological probes is that 
they can be found at very high redshifts l Lamb & Reichart 2000 1. 
One may then use high-redshift GRBs data to probe the cosmol- 
ogy. This is performed by using the (r; a9 , L i30 ) and (V, L iao ) re- 
lations to estimate the luminosity distance at higher redshifts. A 
possible method would consist in using GRBs with z < 1.5 to- 
gether with the luminosity distance estimates from supernovae to 
calibrate the luminosity estimator, and then use this calibration to 
find the luminosity distance of GRBs with z > 1.5. It should be 
noted that this method aims to study the luminosity distance at the 
range 1.5 < z < 5. Later it will be shown that adding information 
obtained in the range z < 1.5 is actually crucial to the study of 
dark energy models. 

With this information, one may estimate the luminosity dis- 
tance at high redshifts, and place constraints on the cosmological 
parameters via a standard \ 2 minimization procedure. 

The aim of this paper is to use this method to constrain the 
GCG unification model of dark energy and dark matter. Even 
though it is shown that the optimum redshift range for studying 
dark energy is around z < 2 I Huterer & Turner 2000 1, as already 
remarked, since the GCG also describes dark matter, a higher red- 
shift range might be relevant for a better understanding of the 
model. 

Our study will be performed in three steps. First we build 
a realistic mock distribution of GRBs in redshift and isotropic 
luminosity space (section 2.2). Second, we test the calibration 
procedure to find what improvements might be achieved in the 
future. To do this we shall consider a fiducial set of calibra- 
tion parameters to generate a mock set of time lags and vari- 
abilities for each GRB. We shall then perform a \ 2 fit to this 
mock data to study the calibration precision. The last step con- 
sists in employing this method to probe GCG models. This 
will be done in a fashion similar to what was already per- 
formed with SNe la iGoliath et al. 200TllWeller & Albrecht 20021 
IDi Pietro & Claeskens 20031 ISilva & Bertolami 20031 A fiducial 
cosmological model is assumed, and regions of constant \ 2 w iH 
be plotted around the fiducial set of parameters. 



2.2 Generating a GRB mock population. 

We describe here the process of population generation. We simu- 
late several data sets, which differ only in size. We consider three 
sample sizes. The first sample is composed of 90 GRBs, and is con- 
sistent with the expected number of GRBs with measured redshift 
z > 1.5 that is what the SWIFT satellite is expected to detect in its 
three year mission I Schaefer 2002 1. The second and third samples 
are larger and contain 500 and 1000 GRBs, respectively, and serve 
as best case scenarios to test whether further data might improve 
the results. 

The GRBs in each sample are distributed in redshift and 
luminosity according to the GRB rate history and luminos- 
ity function based on the model of the star formation 2 from 
|Porciani "&* Madau (20021. The energy spectrum of the GRB is as- 
sumed to be a single power law with index -2.5. The flux limit 
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Figure 1. The redshift distribution of a GRB sample population. 



Parameter 


Value 


Bv 


1Q55.32 


Pv 


1.57 


<J V 


0.20 


B T 


10 50 ' 03 




-1.27 


ay 


0.35 



Table 1. Values used in the calibration test. 

of the SWIFT satellite, / > 0.04 photons cm~ 2 s -1 is applied 
to check whether the GRB can be detected. The observed mag- 
nitude is calculated assuming a flat ACDM cosmological model, 
with £3a = 0.7 and Ho = 70 km s" 1 Mpc" 1 . 

We show an example of a generated sample in FigureQ 

2.3 Calibration procedure. 

|Reichart et al. (200T) proposed a relation between the variability 
and isotropic equivalent luminosity of a GRB such that 



Liso — B V V v , 



(4) 



while |Norris et al. (200 ( proposed a similar relation between r; ag 
and Li SO , 



Bv 
lag' 



(5) 



As already mentioned, the first step in testing the calibration 
procedure consists in establishing a fiducial model. |Schaefer ( 2002 1 
used the nine GRBs with available redshifts at the time to calibrate 
these relations, yielding: 

0.03 



1.57 ; 
-1.27 : 



B v = 

B r 



10" 
= 10 5 



(6) 
(7) 



We shall assume these values as our fiducial model, i.e. we suppose 
that the calibration relations are faithful, and that they are described 
by this set of parameters. For each GRB of the mock luminosity 
distribution, generated as we explained above, we compute the cor- 
responding time lag, r; a9 , and variability, V, through 

1 



log Ti a g = log B T + log L + Random(a T ) 

Pr 

log V — log B v + -j- log L + Random(a v ) 

Pv 



(8) 
(9) 



where the Random(a) term is a pseudo random number drawn 
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Ngrb used in calibration 




a 0v 


o>(mag.) 


40 


0.053 


0.021 


0.68 


100 


0.033 


0.014 


0.66 


200 


0.023 


0.010 


0.66 



Table 2. Results from calibration test. The last column refers to the uncer- 
tainty of the distance modulus determination using our method. 



from a normal distribution with zero mean value and variance a. 
That is, we assume that Ti ag and V have lognormal error distribu- 
tions, with variance a T and a v , respectively. 

The values we use for a v and a T are based on Schaefer (2002 1 
and are essentially dominated by an intrinsic (statistical) error. The 
used fiducial values are shown in Table 1 . 

In the real situation, we would assume that the luminosity dis- 
tance at redshifts smaller than z = 1.5 has been measured indepen- 
dently by SNe la experiments, such as SNAP. Knowledge of the 
peak flux observed on Earth, and the spectral shape N(E) allow 
obtaining the corresponding equivalent isotropic luminosity from 
Eq.0. 

We do not consider errors in this estimate of Li SO , but it should 
be noted that it is expected that the luminosity distance uncertainty 
will be around <j\ og d L ~ 0.01 for future experiments such as SNAP 
<Goliath et al. 2001llWeiler & Albrecht 20021 This translates into 
o"io g L ~ 0.02, a much smaller error than the intrinsic scatter of 
GRBs. Thus, this source of uncertainty should not influence our 
conclusions. 

Notice that we have used low-redshift GRBs from the gen- 
erated sample to calibrate the relations Eqs. and J5) through a 
standard \ 2 fitting procedure. That is, our procedure consisted in: 

• Generating a GRB sample. 

• For each GRB in the sample, assigning a time lag and vari- 
ability via Eqs. {8) and j9j- 

• Fitting the calibration relations Eqs. and ^5) for the low-z 
GRBs from the mock sample. 

• Repeating the steps for different size samples, and study the 
precision of the attained fits. 

Our results are shown in Table 2. The smaller sample size cor- 
responds to what is expected to be found by the SWIFT satellite in 
its three year campaign, the second sample size corresponds to the 
number of low-z GRBs that should be found if there where about 
500 GRBs in total, while the third sample corresponds to about 
1000 GRBs in total. 



2.4 Estimating L iso and the x 2 test. 

For the remaining GRBs that were not used in the calibration pro- 
cedure, the calibrated relations will be used to find Li SO from the 
values of r and V. Thus we obtain two estimates: 



log L v = log B v + p v log V, 

log L T = log B r + /3 r log T iag . 



(10) 

(11) 



These two estimates are combined as weighted averages to produce 
one estimate of Li so , 



log L iao = - ( log L v H — log L T 

2 \Vlo g L v °"logL T 

The next step consists in defining the \ 2 function, 



(12) 



X (p) = 



Ngrb 

E 



Mi — 25 — 5 log d L (zi,p) 



(13) 



where rm is the observed magnitude, A'U is the absolute magnitude 
estimated from Eqs. llOi and and d,L(zi,p) is the standard 
luminosity distance as function of redshift z and the cosmolog- 
ical parameters p. The denominator <r M is the uncertainty in the 
determination of the distance modulus, m* — Mi, This uncertainty 
was calculated using Gaussian error propagation, and assuming that 
there were no correlation terms to be taken into account. The values 
of <7 M we have found are shown in Table 2. We have to minimize 
this function and draw the x 2 contours in order to find the confi- 
dence regions. 

In here we follow the method developed for the study 
of SNe la tGoliath et al."200Tl IWeller & Albrecht 20021 
IDi Pietro & Claeskens 20031 ISilva & Bertolami 2003L We 
use the log-likelihood function \ 2 to build confidence regions 
in the parameter space. In order to perform this one chooses a 
fiducial model, denoted by the parameter vector p/id, and then the 
log-likelihood functions x 2 are calculated based on hypothetical 
magnitude measurements at the various redshifts. The \ 2 function 
is then given by 



N G Rl 



x 2( p j _ f 5logd L (zi,Pfi d ) - 5logd L ( Zl ,p) 



(14) 



Of course, there is no need to minimize the \ function since its 
minimum will correspond to p — Pfid, thus we only need to find 
the x 2 contours corresponding to the desired confidence level. 



2.5 The models. 

The GCG model smoothly interpolates between a dark mat- 
ter dominated time in the past, to an accelerated de Sitter 
phase in the future. Thus, this is a setting that on large scales 
agrees with the observed expansion history of the Universe 
tKamenshchik et al. 20011 IBento et al. 2002b> . The GCG density 
may be written as a function of redshift as 



Pah (z) = PeH , [A, + (1 - A s )(l + z) 3 ^] V(1+Q) 
where p c h,o is the present day density of the GCG, and 



A s — Ap chQ , 



(A + B) 1/(1+a) 



Pchfi 



(15) 



(16) 



where B is the integration constant that appears in Eq. J2j. For 
2 > Owe have a matter dominated Universe, 



Pch (z^>0) = (l~A 3 ) 1/{1+a) (l + z) 3 



(17) 



while in the far future, z = — 1 the GCG behaves as a vacuum 
dominated Universe. 

The GCG unifies dark matter and dark energy, but does not 
take into account the presence of baryons and radiation. The bary- 
onic component has to be considered when studying the implica- 
tions of the GCG model with regards to observational tests, but for 
geometric tests, such as the magnitude-redshift relation we con- 
sider here, the effect of these components is negligible. Therefore, 
throughout the paper we disregard the presence of baryons and ra- 
diation. This does not affect any of our conclusions. 

For the purpose of comparison we choose the so-called 
XCDM model, which consists of two components, cold dark mat- 
ter and some form of dark energy which has a constant negative 
equation of state, w = p/p. We shall use this parameterization to 



test GRBs as a probe of an unspecified dark energy component. Al- 
though this parameterization is not suitable for general dark energy 
models, since in most cases the equation of state changes with time, 
it is adequate to test the cosmological constant model, w = — 1, 
against other models. If one finds that w — —1 is disfavoured by 
the data, then there is a strong indication that the dark energy com- 
ponent is more complex than expected. 

For many years cosmologists were mainly concerned with 
models to which —1 < w < —2/3, however, more recently mount- 
ing evidence that w < — 1, the so-called phantom dark energy 
models tCaldwell 20021 |Carrol, Hoffman & Trodden 2003) , is be- 
ing encountered. In many models there are several theoretical rea- 
sons not to consider w < — 1, most notably that this would lead to a 
breakdown of the weak, the dominant and the strong energy condi- 
tions (p > and p + p > 0, p > \p\, and p + Sp > 0, respectively). 
The weak and strong energy conditions are ones to hold in proving 
well known singularity theorems, while the dominant energy con- 
dition guarantees the stability of a component, such that there is no 
creation of energy-moment from nothing (Hawking & Ellis 1973 1. 

Another feature of some phantom models is that they exhibit 
a future blow up of the scale factor in a finite time, often referred 
to as the 'Big-Rip'. Since such a Universe has a finite lifetime, it 
has been argued that these phantom models solve the coincidence 
problem, that is, they explain why dark matter and dark energy den- 
sities are of the same order at the present time. For a ACDM Uni- 
verse, a long matter dominated period is followed by a rather quick 
transition to a phase such as the one the Universe is now, where 
matter and dark energy have approximate densities. This phase is 
then followed by an eternal vacuum dominated exponential expan- 
sion. Within this framework, the probability of finding ourselves in 
this intermediate and temporary epoch is very small. For phantom 
models on the other hand, the accelerated expansion is not eternal. 
Since the Universe has a finite lifetime, the probability of living in 
the epoch of matter-energy approximate equality is larger than that 
of ACDM models, thus, the cosmic coincidence is not as unlikely 
as for ACDM models. However, these models fail to explain why 
the dark energy component did not start to dominate the evolution 
at a earlier time, namely, why dark energy started dominating the 
cosmic evolution only after the large scale structure had time to 
evolve deep into the non-linear regime, even though in the context 
of the CGC this is exactly what is found l Bento et al. 2004 1. 

Furthermore, besides these theoretical features, the latest 
SNe la data does seem to favour a phantom energy compo- 
nent. Thus, in what follows, we shall consider a phantom model 
for comparison with the GCG model. Moreover, as shown in 
IBertolami et al. (2004) and |Bento et al. (2005t , for z < 2, the GCG 
model is degenerate with phantom models with suitable parame- 
ters. 

As mentioned, we consider two different cosmological mod- 
els. The flat GCG model, which unifies dark energy and dark mat- 
ter into a single component, and the XCDM, which parameterizes 
dark energy in terms of a constant equation of state, w — p/p. 
We consider that the Universe is flat in both models. The GCG 
model is then described by the exponent a, and the quantity A s , 
while the XCDM is described by the parameters w and fi m , the 
non-relativistic dark matter density relative to the critical one. We 
considered two fiducial models. Model I assumes that the Universe 
is described by a GCG model with 1 — A 3 — 0.3 and a = 1. Model 
II assumes the corresponding degenerate Universe, described by a 
phantom XCDM model, with w = —1.4 and Q m = 0.45. 
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Figure 2. Effects of the improvement in calibration of GRBs for CGC mod- 
els. The three sets of curves are indistinguishable. 
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Figure 3. Effects of the improvement in calibration of GRBs for XCDM 
models. The three sets of curves are indistinguishable. 



3 RESULTS. 

We first examined the effect that an improvement in the calibration 
might have in constraining the models. To that end we consider that 
the high-redshift GRBs sample is comprised by 90 GRBs with red- 
shifts z > 1.5, and use the values shown in Table 2 to calculate 
the uncertainty in the determination of Li SO , The effect that an im- 
provement in calibration has in the quality of the cosmological tests 
is shown in Figures|2|and|3] 

We then tested the effect of an increase in the number of high- 
redshift GRBs in the test. We assume for the purpose of calibration 
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Figure 4. Effects of the increase in the number of GRBs for CGC models. 
The curves correspond, from the outer to the inner ones, to 90, 500 and 1000 
high-redshift GRBs 



Figure 6. Effect of adding 100 low-redshift GRBs to a sample of 500 high- 
redshift GRBs for CGC models. Both sets of curves are indistinguishable. 




Figure 5. Effects of the increase in the number of GRBs for XCDM models. 
The curves correspond, from the outer to the inner curves, to 90, 500 and 
1000 high-redshift GRBs. 
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Figure 7. Effect of adding 100 low-redshift GRBs to a sample of 500 high- 
redshift GRBs for XCDM models. The outer curve corresponds to 500 high- 
redshift GRBs, while the inner one corresponds to the 100+500 GRB case. 



that 100 low-redshift GRBs are known, and use the corresponding 
calibration precision shown in Table 2. We next considered larger 
samples of GRBs. Samples with 90, 500 and 1000 high-redshift 
GRBs are examined, and the respective confidence regions are ex- 
hibited. These results are shown in Figures|4]and[5] 

Finally we study the effect of the distribution of GRBs in red- 
shift space. In order to perform this analysis we added 100 GRBs 
with z < 1.5 to a sample of 400 high-redshift GRBs, and found the 



corresponding confidence regions. The results are shown in Figures 
0and|71 

The first conclusion we can draw is that GRBs are not quite 
suitable as stand-alone probes of dark energy. This is seen from 
the lack of impact that an increase in the redshift range has on the 
results in terms of precision and discriminating power. 

The main source of uncertainty comes from the intrinsic sta- 
tistical scatter of the GRB population. As may be seen from Figures 
|2|and|3| calibrating the (r, Li SO ) and (V, Li so ) relations with more 
than about 30 low-redshift GRBs does not substantially affect the 
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results, however an increase in the size of the higher redshift GRB 
population will improve the quality of the constraints that can be 
imposed. This can be seen in Figures [4] and [5] Using more precise 
luminosity estimators, such as the one we will use in the follow- 
ing section, might improve the outcome, even though better results 
can be achieved by probing a greater sample of z < 1.5 redshift 
sources. This is evident for XCDM models, as shown in FigureQ 
Adding more 100 GRBs with z < 1.5 to a sample of 400 high- 
redshift GRBs is more effective than adding more 500 GRBs with 
z > 1.5. 

Regarding the results for each model, the first conclusion is 
that the XCDM model is better constrained than the CGC model. 
It can be seen that, no constraint in the a parameter for the GCG 
model can be imposed, even though an upper limit for A s can be 
obtained. 

As for XCDM models, we point out that there is some poten- 
tial to probe Q m , but the prospect of using this test to probe the 
nature of dark energy is very limited. A upper value for w can be 
obtained, however no lower limit. It should be noted that other tests 
such as SNe la are capable of imposing tighter constrains on fi m 
than the ones we find. 



4 THE GHIRLANDA RELATION. 
4.1 Description. 

The (t, Li 30 ) and (V,Li ao ) relations are notoriously affected 
by the large intrinsic scatter of the data set, which greatly hin- 
ders their use as precision cosmological probes. Quite recently, 
Ghirl ancTa et al. (2004a) have found a surprisingly tight correlation 
between the peak energy of the 7-ray spectrum, E pea k (in the 
v — vF v plot), and the collimation corrected energy emitted in 7- 
rays, £ 7 , for long GRBs. This collimation corrected energy mea- 
sures the energy release by the GRB taking into account that the 
energy is beamed into a jet with aperture angle 9. 

Let us denote by Ei SO the isotropically equivalent energy, in- 
ferred from the isotropic GRB emission. This source frame "bolo- 
metric" isotropic energy is found integrating the best fit time- 
integrated spectrum N(E) [photons cm -2 keV -1 ] over the energy 
range 1 keV - 10 MeV, 



Eiso 



^d'i(z) 



S-yk erg , 



(18) 



where cLl{z) is the GRB luminosity distance, k refers to the k- 
correction (Bloom, Frail & Sari 2001 1 and 5* 7 is the fluence, 



EN(E)dE , 



(19) 



where E is in keV. 

If the 7-ray emission is collimated into a jet with aperture 9, 
then the true energy emitted is E-, — Ei 30 (l — cos 9). Thus, to 
convert Ei SO into E 7 and vice-versa, one needs to know the an- 
gle 9. Under the simplifying assumption of a constant circum-burst 
density medium of number density n, a fireball emitting a fraction 
rj-, of its kinetic energy in the prompts 7-ray phase would show a 
break in its afterglow light curve when its bulk Lorentz factor T 
becomes of the order of T ~ 1/9, with 9 given by <Sari 1999V 



9 = 0.161 



t, 



1 + 2 



3/8 



Eiso.s: 



1/8 



(20) 



where tj et is the break time in days, and Ei ao j,2 = Ei ao j\Q 



-ji J Percentage of y -g-t J 



(f 2 



10% 

5% 

8% 
5% 
17% 
50% 
20% 



0.0027 

0.0006 
0.0027 
0.0025 

0.0259 

0.0069 

0.0081 



5% 

1% 
5% 
5% 
52% 
14% 
16% 



Table 3. Error budget for the Ghirlanda test. The uncertainty in the circum- 
burst density is assumed to be 50%, while the other values are taken from 
Xu (2005). 



The Ghirlanda relation (Ghirlan da et al. 2004al IXu 2005 \ is 
then expressed as 



En, 



10 50 erg 



C 



UOOkeVy 



(21) 



were a and C are dimensionless parameters. 

Using this result, together with Eqs. Jl 81 . 1201 and the defini- 
tion of _E 7 , one obtains an equation for (1l{z), 



C 



En 



100 keV 



= (!-■ 



/ Ei so \ 

Viosoy 



(22) 



Note that both E iso and 9 depend on d^ix), as may be seen from 
their definitions. Solving this equation for <1l(z) we find an esti- 
mate of the luminosity distance at redshift z in terms of the observ- 



ables E p , tj et , n, and fluence 



The error budget for d,L is then found to be ( Xu 2005 1 



m 1 



VS., 



+ 



( 1R. 

\ C 

Sat, 



+ 



E° bs 



+ 



+ 



4(1- VC~e) 2 



+ 



+ 



( cr a E p 
a — m — — 
V a 100 



7? 7 (1 - 77 7 ) 



V 



+ 



Cg 



4(1- VCe) 



(23) 



where Cf = 9 sin 0/(8 - 8 cos 9). 

The distance modulus is given by ^l oos = 5 log dL/10pc, 
thus its error is <r M = 5/ In 5((Jd-L/di / ). We use the error estimates 
from ( Xu 2005 1, which are consistent with what is found in 
the literature jGhirlanda et al. 2004a1 IGhirlanda et al. 2004bl 
IFriedman & Bloom 20051 IGhisellini et al. 20051 

|Xu, Dai & Liang 2005|IM6rtsel & Sollerman 20051 . 

These values yield an uncertainty in /j, t s of order cr M ~ 0.5, 
which is just slightly smaller than the error bars from using both the 
(Li 30 , V) and (Li ao , tl) relations. As may be observed in Table 3, 
the smaller intrinsic scatter of the Ghirlanda relation is balanced 
by its dependence on poorly constrained quantities. Improving the 
calibration will not solve the problem, since the main sources of 
uncertainty are due to the determination of the peak energy E p , the 
jet break time tj e t and the value of the circum-burst density n. 
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Figure 8. Same as Figure[5] but using the Ghirlanda relation. 



-3 




Figure 9. Same as FigureQ but using the Ghirlanda relation. 



5 DEGENERACY BETWEEN MODELS. 

The CGC and the phantom XCDM models are degenerate 
at redshifts z < 2, as shown in Bertolami et al. (2004J and 
|Bento et al. (2005) . Consider the Taylor expansion of the luminos- 
ity distance as 



d L 



-A 



qo)z 2 - - (l - go - 3<?o + jo) 



1 + (1 



go) 



l 

6 

5 3 
o<7o 



12 



gojo 



ko 
24 
+0{z t 



)}, (24) 

were go is the deceleration parameter, related to the second deriva- 
tive of the expansion factor, jo is the so-called "jerk" or statefinder 
parameter l Alam et al. 2003 Sahni et al. 2003 1, related to the third 
derivative of the expansion factor, and fco is the so-called "kerk" 
parameter, which is related to the fourth derivative of the expan- 
sion parameter ( Visser 2004 Dabrowski & Stachowiak 2004 1, all 
evaluated at present. These quantities are defined as 



1 d 2 a 

1 d 3 a 

J(t) = +- -TT3 

a at 1 * 

1 d 4 a 



Ida]- 2 , 
a dt J ' 

1 da] -3 
a dt . 

1 dffl-|- 4 _ 
a dt . 

and are related to each other through the relations 

3 ( p(z) 



k(z) 



,w +1 

q(z)+2q 2 (z: 

dj 
dz 



dq 



(*); 



(*) - Viz) 



dz 



(25) 



(26) 



(27) 



(28) 



(29) 



(30) 



were p(z) and p(z) refer to the total density and pressure of the 
Universe, respectively. With these expressions one can write the 
present day values of the parameters as a function of the cosmolog- 
ical parameters of the model under examination. 
For the XCDM model, we have 



XCDM 

g D 

dg 

dz 

dz 





XCDM 



= -w{l- 



9 2 

= — — w (2 — 3w)(Qn 



1)0, 



(3D 



(32) 



(33) 



while for the GCG model we find 



4.2 Results. 

In Figures [8] and [5] we see the improvements that can be obtained 
by the use of the Ghirlanda relation. We have assumed two redshift 
distributions; one made up of 500 GRBs with z > 1.5, and the 
other made up of 100 GRBs with z < 1.5 plus 400 GRBs with 
z > 1.5. We find that the conclusions one can draw remain essen- 
tially unchanged. Moreover, we see that the high-redshifts of GRBs 
are not quite suitable to adequately study the GCG model, while it 
allows for the XCDM model to place a rather useful limit to the to- 
tal amount of matter, although not much can be learned with respect 
to the dark energy equation of state. 



GCG 
1o 



= 2< 1 - 



dq 
dz 

dz 



GCG 


GCG 



9 



A s )(l + a), 



l + -a(l 



As ) As 



(34) 



(35) 



(36) 



For the redshift range probed by SNe la one may neglect 
terms beyond the cubic power in redshift in Eq. 1241 . SNe la data 
indicate that these models have the same deceleration and jerk 
<Bertolami et al. 2004IIBento et al. 2005> , that is, they are degener- 
ate for the considered redshift range. Thus, imposing this equality 
one finds the relationship between parameters 
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w = a(A s - 1) - 1, 



(l + g)(l- A a ) 
l + a(l-A s ) ' 



(37) 



(38) 



This degeneracy holds for SNe la, for the maximum probed 
redshift of about z ~ 2. As the redshift range allowed by GRBs 
is greater, one hopes to test higher order terms in Eq. I24i . At next 
order, the degeneracy is broken by the kerk parameter, that is, even 

■r „GCG XCDM -GCG -XCDM „ fi „ . ,, . ,GCG i 

it <7o = <?o . Jo = Jo . one hnds that k 
j^xcdm p Qr mstance; jf one considers a = 1 and j4 s = 0.7, then 
kQ Ca = —0.3, while for the corresponding XCDM model (with 
w = -1.3 and Q m = 0.46), the kerk is ko CDM = -4.27. 

This procedure can be viewed as a consistency test. Indeed, 
consider the GCG model. Once the deceleration, the jerk and the 
kerk of the Universe are measured, one can extract the values of a 
and A s . These values yield a value for the kerk, call it fcg, which 
can be compared to the measured value of the kerk, ft™. Whether 
both values are close or not, one may reach a conclusion regarding 
the suitability of the model. The same reasoning can be applied to 
the XCDM model. This is, of course, a simplification, since one 
has to consider the accuracy of measurements and the statistical 
significance of each of the parameters in Eq. 1241 . as well as to 
consider correlations between them. 



6 DISCUSSION AND CONCLUSIONS. 

We have verified that the medium and high redshift of GRBs do 
not impose too strong constraints on the GCG model. At such high 
redshifts, the expansion factor for the GCG becomes 



H ch (z » 0) = Q ch (l - A s ) 



l/(l+c) 



(39) 



where Q. c h is the CGC density relative to the critical one. Since 
we are considering a flat Universe made up mostly by the GCG, 
£l c h = 1. Thus, the luminosity distance at such high redshifts de- 
pends essentially on (1 — A s ) 1 ^ 1+OL \ It is found that the confi- 
dence regions for the GCG model approximately follow the line 
(1 — A s ) 1 ^ 1+a ' ) = const. Hence, since A s and a are strongly 
correlated by this expression, neither one is constrained. 

As for the XCDM model, a high-redshift population of 
GRBs is suitable to constrain the total amount of matter, Q m , 
but are poor probes of the dark energy equation of state. This 
does not mean that GRBs have no use for the study of dark en- 
ergy, or that expanding the Hubble diagram is meaningless. It 
is a common trend in the recent literature to argue that a pre- 
cise prior estimate of Q, m enhances the ability of SNe la to 
constrain the dark energy equation of state l Goliat het al. 200 Tl 
IWeller & Albrecht 2002llDi Pietro & Claeskens 20031 . This is also 
true for GRBs. A high-redshift sample will constrain fi m , while 
a low-redshift sample will place constraints on the equation 
of state. A low-redshift sample alone will fail to constrain 
the equation of state since it will be hindered by large un- 
certainties in Q. m . This is somewhat similar to what happens 
in SNe la tests. A low-redshift sample will constrain the nui- 
sance parameter A4, while the high-redshift sample constrains 
the cosmological parameters jP admanabhan & Choudhur y 2003 1 
|Choudhury & P admanabhan 2005)] 

Our study also reveals that the effectiveness of the GRB test 
depends on the model under examination. GRBs allow imposing 
some constraints on the XCDM model, but very little can be ex- 
tracted for the GCG model. This may be explained by two factors. 
The first was explained above, and is related to the fact that the 




Figure 10. Total equation of state of the Universe as a function of redshift. 
The value z = — 1 corresponds to the very far future. The dashed line 
represents our GCG fiducial model, while the solid line corresponds to the 
XCDM phantom model. 



region probed by GRBs depends on an unique function of both pa- 
rameters. The second is related to how late the Universe stops being 
matter dominated. As shown in Fieure flOl the CGC model stays in 
the matter dominated phase for a longer time, thus the period of 
time during which the dark energy behaviour affects the cosmic 
expansion is slightly larger for the XCDM model. Since GRBs es- 
sentially probe larger redshifts, models in which the transition be- 
tween the matter dominated and accelerated expansion is shorter 
cannot be well constrained by this test. By the same reasoning, 
GRBs might be actually useful in testing models where the tran- 
sition between matter dominated expansion and dark energy driven 
acceleration started at larger redshifts. 

Thus, we find that despite all intrinsic limitations, larger sam- 
ples of GRBs can potentially determine the total amount of matter 
in a XCDM model. This independent determination of il m will be 
useful in SNe la studies of the dark energy equation of state. 

Furthermore, GRBs may also be used to break the degeneracy 
between models, since at high-redshifts the effect of higher order 
terms in the Taylor expansion of dh (z) must be taken into account. 
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